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ABSTRACT 

We use quantitative metrics to characterize the variation of C IV A 1549 broad absorption lines 
(BALs) over 3-6 (rest-frame) years in a sample of 13 quasars at 1.7 < z < 2.8 and compare the results 
to previous studies of BAL variability on shorter time scales. The strong BALs in our study change in 
complex ways over 3-6 yr. Variation occurs in discrete regions which are only a few thousand km s 
wide, and the distribution of the change in absorption equivalent width broadens over time. We 
constrain the typical C IV BAL lifetime to be at least a few decades. While we do not find evidence 
to support a scenario in which the variation is primarily driven by photoionization on multi-year time 
scales, there is some indication that the variation is produced by changes in outflow geometry. We do 
not observe significant changes in the BAL onset velocity, indicating that the absorber is either far 
from the source or is being continually replenished and is azimuthally symmetric. 

It is not possible in a human lifetime to expand the time scales in our study by more than a factor 
of a few using optical spectroscopy. However, the strong variation we have observed in some BALs 
indicates that future studies of large numbers of BAL QSOs will be valuable to constrain BAL lifetimes 
and the physics of variation. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: individual (LBQS 0010-0012, LBQS 
0018+0047, LBQS 0021-0213, LBQS 0051-0019, LBQS 0055+0025, LBQS 0109- 
0128, LBQS 1208+1535, LBQS 1213+0922, LBQS 1234+0122, LBQS 1235+1453, 
LBQS 1243+0121, LBQS 1314+0116, LBQS 1331-0108) — quasars: absorption 
lines 



1. INTRODUCTION 

Broad Absorption Line Quasars (BAL QSOs) exhibit 
broad, ultraviolet (UV) line absorption troughs spanning 
thousands, or tens of thousands, of km s _1 . By studying 
the evolution of these dramatic features, we hope to learn 
about the structure and dynamics of absorbing material 
in active galactic nuclei (AGNs). Because QSO redshifts 
of z ~ 2 are needed to shift high-ionization BAL lines 
into the optical waveband, observations of BAL QSOs 
must span decades in order to characterize multi-year 
BAL variability for high-ionization lines such as C IV 
A1549. 

BAL troughs are p resent in about 20% of QSOs (e.g., 
iHewett fc Foitzll2~003D . produced by absorption from lines 
such as C IV A1549 and Mg II A2799. BAL absorption 
can extend to very high velocities, although above ks 0.1c 
absorption from C IV A1549 becomes confused with lines 
from Si IV and O IV. BAL QSOs with broad absorption 
from low ionization stages such as Mg II or Al III are 
classified as "low-ionization BAL" (LoBAL) QSOs, while 
those with absorption only from ions at higher ionization 
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stages are called "high-ionization BAL" (HiBAL) QSOs. 

It is commonly thought that BAL outflows are seen 
when QSOs are observed at large inclination angles with 
respect to accretion disk. In this model, the line of sight 
to the central source passes through an equatorial disk 
wind, which imprints a broad absorption signature on 
the continuum. BAL outflows would be present in ev- 
ery QSO, but only cover about 20% of the QSO sky. 
Modeling this scenario has successful ly reproduced many 
properties of BAL absorption (e.g.. iMurrav et all 119951 : 
iProga et alJl200C( ). Other studies have modeled the BAL 
outflow as an orientation-independent evolutionary ef- 
fect, caused b y the expu lsion of a thick shroud of gas 
and dust (e.g.. lVoit et~al]ll993t [Gregg et al.ll2006f ) which 
may have been deposited in the nucleus by galaxy colli- 
sions (e.g., Canaliz o fc Stocktonll200lHJunkkarinen et all 
l200ll lGregglt^[200l . 

BAL outflows have a complex structure. The covering 
factor, or geometric fraction of the emission which is ob- 
scured by the absorber(s), can vary strongly across the 
absorption trough a s a function of outflow velocity (e.g., 
Ide Kool et alJ feOOl). As a result, the absorption is often 
saturated but not black, so that the outflow geometry 
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(with respe ct to the emitter) determines the absorption 
profile (e.g.. lArav et al.lll999al ). The ionization level also 
apparently varies with velocity, and, at least in the case of 
PG 0946+301, the covering fa ctor increases with the ion- 
ization level of the absorber (jArav et~allll999bh . Spec- 
tropolarimetric studies have found that the emission at 
the bottom of deep troughs is highl y polarized and li kely 
scattered around the absorber (e.g. lQgle et al.llT999D . 

In this study, we take advantage of the ov erlap between 
the L arge Bright Quasar Survey (LBQS; iHewett et al.l 
119951 and refe rences therein) an d the Sloan Digital Sky 
burvey (SDSS; lYork et al.ll2000h to examine variation in 
a sample of 13 BAL QSOs over rest-frame time scales 
of 3-6 yr. 1 Previous studies have focused on individ- 
ual sources (j jHH or on a larger number of sources over 
shorter (< 1 yr) time scales f m.3p . For comparison, we 
plot BAL equivalent widths (EWs) and times between 
observa t ions i n Figure [1] for our st udy and the samples of 
iBarlowi (|1993l hereafter B93) and lLundgren et all (|2007l 
hereafter L07). 2 In this paper, we adopt the conven- 
tion that absorption features have negative EWs, and 
the EWs refer to rest-frame values. 

BAL variation may be driven by several processes 
working on different time scales. For sufficiently dense 
absorbers, photoionization time scales can be weeks or 
less (e.g., iKrolik fc Krissl l2001h . Flow time scales for 
structural variation can be much longer, and proceed in- 
dependently of the continuum. For example, at a speed 
of 10,000 km s _1 , it would take «3 yr to cross a re- 
gion of size ~ 10 17 cm (~ 10, 000 gravitational radii for a 
10 s M Q black hole). Eventually, surveys of large numbers 
of BALs should be able to distinguish variation trends at 
different time scales. This study takes one step toward 
that goal by characterizing multi-year variation in a sam- 
ple of C IV BALs in the context of previous studies on 
shorter time scales. 

In the remainder of this introduction, we briefly re- 
view studies of BAL variation for individual BAL sources 
( fll-l[) ; for narrower, "mini-BAL" absorbers (" £1 1 . 2[> ; and 
also for samples of objects with BAL absorbers ( §1.3p . 
Throughout this paper, we use a cosmology with Hq = 
72, Q M = 0.3, and Q A = 0.7. 

1.1. BAL Variation in Individual Objects 

In this section, we summarize recent observations 
of BAL variability in individual sources to give con- 
text for our study. High-velocity absorption compo- 
nents have va ried over a few months to a few years in 
some sources dFoltz et al.lll987btlSmith fc PenstorJ [T98l 
iTurnshek et al.lll988HVilkoviskii fc Irwinl l2001) . In other 
cases, the variation has occurred acro ss a wide range o f 
velocities and ionization stages (e.g., lArav et "all l2001h 
or in multiple absorption troughs (e.g. jHall et al.l l2002). 
Dramatic X-ray a bsorption variation has been observed 
in PG 1004+130 (iMiller et al.l l2006l) and PG 2112+059 
(jGallagher et al.ll2004D . while UV absorption also varied 
subtly in the latter source. 

Some QSOs have varied so strongly that the BAL sta- 
tus of the object has changed. The Mg II absorber of 

1 Throughout this work, all durations are given in the rest frame 
defined by the QSO redshift, unless otherwise noted. 

2 The method of calculating EWs varies slightly among the stud- 
ies, but the differences should not strongly affect comparisons be- 
tween the samples. 



LBQS 0103-2753A disappeared over 6 yr, changing the 
techni cal classification of that source from LoBAL to Hi- 
BAL (punkkarinen et al.|[200ll) . Narrow C IV and Si IV 
absorption lines in the radio-loud QSO TEX 1726+344 
broadened over 3.5 yr, cha nging its classification from a 
non-BAL to a BAL QSO (|Mall2002f ). 

Photoionizat ion has been propo sed as a driver of BAL 
variation (e.g.. iBarlow et al"Ill989l ). but this explanation 
can be problematic. In several cases, abso rber variation 
does not appear to track the c ontinuum (|Barlow et al.l 
Il992t iMichalitsianos et al.|[l996[ ). The absorber may be 
too tenuous to change its ionization state on continuum 
variation time scales, or BAL variation may be governed 
by additional structural factors. 

Despite the high velocities and complex structures 
commonly observed in BALs, there have only been 
a few cl aims for detections of ac celeration in BAL 
troughs. IVilkoviskii fc Irwin] ((2001) presented a spec- 
trum of Q1303+308 which showed strong BAL varia- 
tion over 6 yr, particularly in the highest-velocity com- 
ponents. They also reported that the outflow velocity 
had increased in magnitude by 55 km s -1 , correspond- 
ing to an acceleration of «0.03 cm s~ 2 . Rupk eet al.l 
(2002) point out that one high-velocity component of 
the Na I D line of the near by QSO Mrk 231 accelerated 
over about 10 yr. Recently, Hall et al. (2007) compared 
SDSS and ESO Very Large Telescope (VLT) spectra of 
SDSS J024221.87+004912.6, finding evidence for an av- 
erage acceleration of the C IV BAL trough of 0.15 cm s~ 2 
over 1.4 yr. 

1.2. Variation in Narrower Absorption Systems 

Narrower intrinsic absorption lines share many p roper- 
ties in common with BALs (e.g.. lArav et al.lll9 99b). so it 
is worthwhile also to consider the variation of mini-BALs 
(i.e., broad absorption troughs which are somewhat nar- 
rower than the 2000 km s _1 required for a BAL). Mrk 231 
shows broad (« 1500 k m s -1 wide) absorption lines from 
several different ions (Bokscnb erg et ail 119771 ) and has 
been the subject of many investigations. A new velocity 
component for Na I an d He I appeared betw een 1984 De- 
cember and 1988 May (|Boroson et alJI 19911 ). Na I D line 
components have continued to vary over the last decade, 
and the optical polariz ation of this source h as also been 
observed to vary (e.g..lGallagher et alJ [2005h. 

Recently, iMisawa et al.l ((20071 ) reported monitoring of 
the QSO HS 1603+3820 (z = 2.54) which hosts a mini- 
BAL. They tracked substructure variation in the mini- 
BAL trough over about 1.2 yr at high spectral resolution 
with the Subaru High Dispersion Spectrograph (HDS). 
The EWs of individual components of the mini-BAL 
varied together, and this variation was correlated with 
changes in the covering factor of the absorber. 

1.3. Studies of BAL Variability in Samples 

In an early study of BAL variability, B93 monitored 
spectra of 23 BAL QSOs over times up to 1 yr. BALs in 
four QSOs changed dramatically during the study, while 
those in 11 more objects changed at least marginally. 
The BAL features varied in intensity rather than in ve- 
locity (within the limits of the data), and there was some 
tendency for BAL changes to correlate with broad-band 
continuum variation. These findings supported the hy- 
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pothesis that photoionization is the primary driver of 
BAL variability. 

Recently, L07 searched for variability in the spectra of 
29 BAL QSOs which had been observed by the SDSS in 
two epochs separated by up to 4 months. They allowed a 
somewhat looser definition of "BAL," retaining in their 
sample objects with absorption velocity widths down to 
1000 km s" 1 . They found that the magnitude of the frac- 
tional change in EW, \AEW/ (EW)\, was greatest over 
the longest times between epochs and for the weakest 
BALs. The strongest fractional changes in EW occurred 
at outflow velocities of -12,000 to -16,000 km s _1 . While 
half of the varying absorption troughs showed changes in 
velocity width, variation in trough depth was more com- 
mon. No correlation was found between changes in the 
BAL EW and UV continuum variability. 

2. BAL SAMPLE SELECTION 

The LBQS contains spectra of 1055 QSOs over a red- 
shift range 0.2 < z < . 3.4 obtained between 1986 and 1989 
(|Hewett et al.lll995l ). The resolution of the spectra used 
in the current study is «6 A (FWHM), with a sampling 
of «2.5 A per bin (in the observed frame) (|Morris et al.l 
1991). The LBQS spectra were kindly provided by P. 
Hewett for this st udy. BAL QSO s in th e LBQS have 
been cataloged by iHewett fc Foltd (|2003h : we refer to 
these objects as "HF03 BAL QSOs." 

The SDSS QSO spectra have been obtained since the 
year 2000 at a spectral resolution of «3 A with wl A 
per spectral bin (in the observed frame). The SDSS 
QSO catalog for Data Release 5 (DR5) identifies 77,429 
QSOs in SDSS DR 5 and gives spectr oscopic redshifts for 
these sources (|Schneider et alj |2007). We have searched 
in the public SDSS database for objects that correspond 
to LBQS QSOs in order to obtain the maximal set of 
sources which have both LBQS and SDSS spectroscopy 
available. 

Most LBQS spectra were obtained using the Multi- 
ple Mirror Telescope (MMT) spectrograph, with spec- 
tral coverage from 3200 < A < 6400 A. The LBQS 
spectra extend out to « 7500 A, but are contaminated 
by the second-order spectrum for A > 6500 A, and we 
discard this wavelength range. The SDSS spectra cover 
longer wavelengths, 3800 < A < 9200 A. QSO redshifts 
of 1.6 < z < 3.1 are required to fit continua and cover the 
C IV absorption region from 1400-1550 A in both spec- 
tra. We converted the LBQS air wavelengths to vacuum 
wave lengths for comparison with SDSS spectra (|Mortonl 
fl99l . 

There are 127 QSOs (with or without BALs) which 
have full C IV region coverage in both surveys. No ob- 
jects have coverage of both the C IV and Mg II regions 
in both epochs, so we are unable to study C IV BALs 
in relation to Mg II absorption (or lack of absorption). 
We have considered analyzing Mg II BALs individually, 
but several factors complicate tests for broad Mg II ab- 
sorption in our spectra. Broad emission features in the 
Mg II region obscure absorption and the underlying con- 
tinuum, and the spectra do not cover the red side of the 
emission region in many cases. 

We have included one BAL QSO, LBQS 0021-0213, 
which we previously observed with the Hobby-Eberly 
Telescope (HET) Low-Resolution Spectrograph (LRS) 



on 2001 December 10 at a sufficient resolution to in- 
clude in our study. (This source was not observed by 
the SDSS.) The source was observed with a 2" slit using 
the g2 grism and the GG385 filter, giving a resolution of 
8 A. We reduced the spectrum using standard procedures 
with the Image Reduction and Analysis Facility (IRAF) 
version 2.12. It was flux-calibrated using observations of 
a standard star taken the same night. While the abso- 
lute flux measurement for this source is not known with 
sufficient accuracy to be used in our study, the structure 
of broad absorption features (with a continuum divided 
out) is reliably represented in the spectrum. We con- 
verted air wavelengths to vacuum wavelengths as we did 
for LBQS spectra. 

The error i n the LBQS flux c alibration is estimated 
to be «15% (|Hewett et alJll991h . We re-normalize the 
SDSS spectra so that the synthetic filter fluxes match the 
photometric fluxes. We divide out the continuum in most 
of our study, rendering the absolute flux normalization 
unimportant. In cases where we do use the continuum 
flux, we assume the SDSS error is small compared to the 
LBQS error. 

2.1. Continuum Fitting 

In this section, we describe our method of fitting con- 
tinua and emission line profiles to LBQS, SDSS, and 
HET spectra of all 128 sources with two epochs of com- 
plete coverage in the C IV region. We correct all spec- 
tr a for Galactic extinc tion using the reddening curve 
of ICardelli et a l. (1989) with the near-UV extension of 
lO'Donnelll (fl99l ). We obtain E(B-V) from the NASA 
Extraga lactic Database (NE D) , which uses the dust 
maps of ISchleeel et al.l (|1998l ). 

We convert all spectra to the rest frame using the 
SDSS DR3 QSO catalog redshift when available, oth- 
erwise we use the redshift assigned by the SDSS data- 
processing pipeline. 4 (These redshifts also agree with the 
redshifts assigned by the LBQS.) Some of the emission 
lines used by the SDSS pipeline to determine redshifts 
may be distorted by BAL features, but it does not ap- 
pear that this effect has significantly impacted SDSS red- 
shift determination for our sources. The emission from 
C III] A1909 is only weakly absorbed, and our fits of this 
line in the SDSS spectra are not greatly influenced by 
Al III A1857 emission 40 A to the blue side. Compar- 
ing the redshift of the C III] emission line in our fits to 
the SDSS systemic redshifts, we find a mean difference 
(Az) = (zsdss — z ciii}) = 0.007 with a maximum differ- 
ence max(|Az|) = 0.014. An RMS difference of 0.006 in 
redshift determination has been obs erved in sources ob- 
served multiple times in the SDSS (jWilhite et all 120051 : 
ISchneider etafll2007f ). Therefore, we do not believe that 
BAL absorption has led to large errors in redshift for our 
sample. 

We fit a continuum to the following regions (when 
present), which are expected to be relatively free of line 
emission and absorption features: 1250-1350, 1600-1800, 
1950-2050, 2150-2250, and 2950-3700 A (rest frame). At 

3 http:/ /nedwww. ipac.caltech.edu/ 

4 The DR5 QSO catalog JSchneider et al.ll2007h had not been 
published when we began our study, but the only source with a red- 
shift different from that of the DR5 catalog is 002444.11+003221.4 
(at z = 0.40) which differs by [Az| = 0.007. This source does not 
affect our study. 
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each iteration, we ignore any spectral bins which deviate 
by more than 3<r from the previous iteration's model fit. 
This allows us to exclude additional regions which con- 
tain unpredictable broad emission or absorption struc- 
tures such as those seen in some cases on the red side 
of C IV A1549 and the blue side of Mg II A2799. We 
fit a power law continuum intrinsically red dened usin g 
the S mall Magellanic Clou d (SMC) curve of HH (fl99l . 
As in iHopkins et all (|2004D . we find that SMC-like red- 
dening reproduces the observed spectra more accurately 
than a Galactic reddening model does. Although the un- 
derlying continuum structure may be more complex than 
a reddened power law (due, e.g., to breaks in the contin- 
uum) , our experience with the sources in this paper (and 
also a large number of SDSS spectra not used in this 
study) indicate that a reddened power law reproduces 
the continuum structure well with a small number of fit 
parameters. 

We use our best-fit continuum model to calculate 
-Fi,(2500 A), the continuum flux density at (rest-frame) 
2500 A. Because broad emission line features can con- 
taminate the spectrum at 2500 A, we prefer to use the 
continuum fit rather than a direct measurement of the 
continuum at that wavelength. This method also allows 
extrapolation of the continuum to 2500 A when that re- 
gion is not present in the observed spectrum. We sim- 
ilarly calculate F„(1400 A), the continuum flux density 
at (rest-frame) 1400 A. For the BAL QSOs in our study, 
the continuum is reasonably well-constrained at 1400 A. 

With the continuum fixed, we fit Voigt emission line 
profiles, corresponding to Si IV A1399, C IV A1549, C III] 
A1909, and Mg II A2799, when these wavelengths are 
present in our spectrum. We use a single Voigt profile 
for each emission line, as additional emission structure is 
not generally required by our data. We fit each emission 
line iteratively. At each iteration we ignore "absorbed" 
bins which are more than 2.5ct below the previous model 
fit. We do not consider the Voigt profiles to be physical 
models. They represent the overall emission line profile 
well with a small number of fit parameters. 

Visual inspection indicates that our fitting algorithm 
generally does an excellent job of reconstructing the un- 
derlying continuum and emission line profiles. In a small 
number of cases, the continuum fits are unsatisfactory. 
We manually adjust the continuum in these cases and 
then re-fit the emission lines as before to obtain a vi- 
sually satisfactory fit. For a few BAL QSOs where the 
SDSS spectrum does not extend blueward of the C IV 
absorption region, we have used the LBQS spectrum as 
a guide to estimate the SDSS continuum. The final con- 
tinuum fits are shown in Figure [2] 

Although we fit a reddened power law model to the 
continuum, we do not interpret our fit results as phys- 
ically meaningful models. BAL QSOs are more intrin- 
sically reddened than non-BAL QSOs (|Wevmann et al.l 
I199H iBrotherton et alj 120011: iReichard et al.ll2003[) . and 
it can be difficult to determine the extent of continuum 
reddening when the shorter-wavelength spectral regions 
are strongly affected by both broad line emission and 
absorption features. Differences in the telescopes and 
calibration used for the LBQS and SDSS further com- 
plicate comparison s between the spectra (e.g., §2.3 of 
iHewett et al1ll991l ). Because the LBQS and SDSS spec- 



tra span different wavelength ranges, we also expect that 
the distribution of fit parameters may differ. The SDSS 
spectra cover a broader range of longer wavelengths, pro- 
viding tighter constraints on the power law continuum 
and weaker constraints on the intrinsic reddening. 

Throughout this study, we refer to the ratio spectrum, 
i?(A), which is constructed by dividing the observed spec- 
tral flux density by the continuum and emission line fit 
model. For the C IV region (1400-1550 A), this includes 
the C IV and Si IV emission lines. In velocity space, 
we call the ratio spectrum R(v). We consider outflow 
velocities and absorption equivalent widths to be neg- 
ative throughout. Changes in measured quantities are 
calculated by subtracting the value in the earlier epoch 
from that of the later epoch. With these conventions, 
AEW < corresponds to an increase in absorption 
strength over time. 

2.2. Identifying BAL QSOs 

BALs are traditionally identified according to the bal- 
nicity index, BI, which approximately expresses the ab- 
sorption EW in km s _1 of a trough which spans at least 
2000 km s _1 and is offset by a " detachment velocity" o f 
at least -3000 km s _1 from rest (Wcvma nn et al.lll991h . 
If no trough meets these criteria, the BI = 0. BAL QSOs 
are defined as those QSOs which have BI > for some 
absorption line. 

The constraint of a detachment velocity 
<-3000 km s _1 may cause one to overlook some 
objects with broad absorption in our sample. We 
investigate this possibility by defining BIq, which is 
calculated similarly as BI but integrated to a detach- 
ment velocity of km s" 1 . We find that we are not 
missing BAL QSOs, as BI > whenever BIq > 
for all the sources. Some BALs do extend beyond a 
velocity of -3000 km s _1 , and BIq > BI in these cases. 
Further modifi cations to the trad itional BI index were 
considered bv lTrump et alJ (|2006h . but BI and BI are 
sufficient for our purpose of simply identifying BALs. 

Our sample of 128 QSOs includes 9 sources which 
were identified as HF03 BAL QSOs. Our continuum fits 
give BI > for 8 of these sources (0018+0047, 0051- 
0019, 0109-0128, 1208+1535, 1235+1453, 1243+0121, 
1314+0116, and 1331-0108). We also include 
1234+0122, which was identified as an HF03 HiBAL 
with a low BI = 4. Our continuum fit technically 
gives BI = for 1234+0122, but multiple strong ab- 
sorption features are clearly present, with one just un- 
der 2000 km s _1 wide. (The absorption appears to 
have weakened by the SDSS epoch.) For three addi- 
tional sources (0010-0012, 0055+0025, and 1213+0922), 
we measure BI > 0, although the absorption is weak in 
these sources and they were not classified as BAL QSOs 
by HF03. The final source, 0021-0213, was observed in 
the LBQS and also with the HET, but not in the SDSS. 
Including 1234+0122, we have 13 QSOs in our sample 
with C IV BALs, i.e., BI > 0. The ratio spectra, R(v), 
for the BAL QSOs in our sample are shown in Figure El 
In this figure, the SDSS spectra have been convolved with 
a Gaussian of the appropriate width to match the reso- 
lution of the LBQS spectra. 

We have assumed that the continuum in the 
1400-1550 A region is a smooth, reddened power law. 
In fact, that may not be the case. While the narrow 
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absorption features in our sample spectra are reasonably 
attributed to absorption, we cannot exclude the possibil- 
ity that broad, shallow features are due to structure in 
the continuum emission. The most plausible candidate 
in our sample for this effect is 1213+0922, which has a 
broad (15,000 km s _1 wide), shallow, wedge-shaped fea- 
ture with a minimum at about -16,000 km s _1 . This 
feature is seen in both the LBQS and SDSS epochs. We 
cannot rule out the possibility that the feature is due to 
continuum structure, but if it is, the results of this study 
would not be significantly impacted. 

None of the BAL QSOs in our sample was considered 
by L07. One source in our sample, 1331-0108, was ob- 
served twice in 1991 by B93. We do not directly compare 
our measurements to those of B93 for this object because 
we cannot account for differences in analysis such as con- 
tinuum placement. We include 1331-0108 in both sam- 
ples when comparing results between studies as if it were 
an independent source in each study. 

The dates of the LBQS and SDSS (or HET) observa- 
tions are given in Table [TJ The measured properties of 
our sources are listed in Tableland Table [3] We discuss 
our BAL QSO sources in detail in the Appendix. 

3. BAL VARIABILITY METRICS 

In this section, we apply various metrics to characterize 
and constrain BAL variation quantitatively. 

3.1. Variation in Velocity Space 

Figure [3] shows the C IV region (1400-1550 A) ratio 
spectra R(v) for the 13 QSOs in our sample. Clearly, 
BAL absorption variation occurs frequently on multi- 
year time scales. In some cases, the entire BAL varies 
(e.g., 0109-0128), while in others, only part of the BAL 
changes (e.g., 0051-0019). In this section, we define a 
criterion to identify velocity regions for which BAL ab- 
sorption has varied between epochs. We consider the 
typical width of a varying region ( £13.1. 1| . the strength of 
variation in these regions ( ^3 . 1 . 2[) . and the dependence of 
variation on the overall absorption strength in a region 
( 35X5)1 . In order to compare the LBQS, SDSS, and HET 
spectra directly, we rebinned all spectra to a common 
grid with bins 1 A wide. This binning oversamples the 
LBQS resolution by ?»2-3 and conveniently allows com- 
parisons between sources and across epochs on a common 
velocity-space grid. For the studies in this section only, 
we smoothed the observed-frame LBQS spectrum with a 
boxcar three bins wide before rebinning in order to re- 
duce scatter without significantly degrading resolution, 
and convolved the SDSS spectrum with a Gaussian in 
order to match the resolution of the LBQS ratio spectra, 
which is «3 A in the rest frame. In this section, we use 
metrics which are relatively insensitive to the measure- 
ment errors in individual bins. 

3.1.1. Characteristic Velocity Width of Variation 

To identify a set of velocity regions in each BAL where 
variation has occurred, we determine the regions in the 
C IV absorption spectrum of each BAL where the ratio 
spectra, R(v), differed by at least 0.1 over 1200 km s _1 
(6 bins) or more between epochs. In bins where the 
ratio spectrum rose above the continuum (R(v) > 1), 
we pegged the value of the ratio spectrum to 1 in order 
to minimize the effects of statistical noise and emission 



structure we have not modeled. While our definition of 
"varying region" allows for cases where AR(v) > 0.1 in 
some bins and AR(v) < 0.1 in other bins of the same 
region, there are only two instances where this actually 
occurs (in 1235+1453 and 1314+0116). This indicates 
that the identified regions are not strongly influenced by 
noise in the spectrum. 

We identify 30 varying regions, and 12 of our 13 sources 
have at least one varying region. The regions are listed 
in Table 3] and are indicated on the ratio spectra in Fig- 
ure [3) In Figure [4) we display the distribution of varying 
regions together with a histogram of the number of times 
a particular velocity was included in a varying region. 
As the figure shows, variation occurs in our sample with 
roughly equal probability from -6000 to -24,000 km s _1 . 

This definition of "varying region" does not identify all 
types of BAL variation. For example, it does not flag the 
trough contraction in 0010-0012, and it is not sensitive 
to acceleration in narrow troughs. It is intended to be 
a well-defined, reliable way to flag significant variation 
without referring to the shapes of individual absorption 
troughs. Averaged across all sources, LBQS error spec- 
trum is about 0.1 in a 1 A bin. Given that the smallest 
varying regions we identify are 6 bins wide, it is unlikely 
that many of them are due to statistical fluctuations. 

Figure [5] shows the number of times variation was ob- 
served in a region with a given velocity width. Varia- 
tion tends to occur on small (<2000 km s _1 ) velocity 
scales. Even the largest variation widths (5000 km s _1 ) 
are narrow compared to the extent of strong BAL ab- 
sorption. The velocity width of a varying region is not 
significantly correlated with the outflow velocity of that 
region, according to a Spearman rank correlation test. 

In a few cases, two varying regions are sufficiently close 
together to raise the possibility that they are separated 
only because statistical noise obscured the variation in 
the intervening bins. Even if these regions were merged, 
it would not affect our conclusion that variation occurs 
on relatively small velocity scales. 

We define Xa« as the sum of velocity widths of all 
varying regions in a single source. For example, if one 
BAL QSO had two varying regions of velocity width 2000 
and 3000 km s _1 , then £ All = 5000 km s" 1 for that QSO. 
The values of Sa„ for our objects are listed in Table 2) 
We do not find a significant correlation between Ea-u and 
At sys , the rest-frame time between epochs, but the range 
of At sys covers only a factor of « 2 in our sample. 

3.1.2. Strength of Variation 

For each velocity region determined to vary ( 33.1.1|) . we 
determine the magnitude of the mean difference, |(Ai?)|, 
between the ratio spectra in that region. (Here AR(v), 
the change between epochs in the bin at velocity v, is av- 
eraged over all velocity bins in a varying region.) Figure[6] 
shows the distribution of |(Ai?)|. Because we imposed a 
threshold of |Ai?| > 0.1 for a bin to be considered part 
of a varying region, we do not expect many cases with 
|(AiZ)| < 0.1. (By our definition, a varying region of 
contiguous bins could, in principle, have |(Ai?)| < 0.1 if 
AR > 0.1 in some bins and AR < —0.1 in others; this 
happens in only two cases.) 

The distribution of |(Ai?)| is strongest at 0.15-0.25, in- 
dicating that changes in absorber depth of |(Ai?)| < 0.25 
are favored on multi-year time scales. We find no signif- 
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icant correlation between | (AR) | and the velocity width 
of the varying region or the average velocity of the vary- 
ing region. (Ail) is positive for 13 of 30 varying regions. 
Surveys with larger samples will be able to constrain fur- 
ther any positive (negative) bias in the distribution of 
(Ail) which could arise if BALs strengthen and weaken 
on different time scales. 

3.1.3. Dependence of Variation on Absorption Strength 

We would like to determine whether the properties of 
variation are dependent on the depth of the absorption 
in the varying region. To do this, we consider three dis- 
tributions. The first, Rail, is defined as the distribution 
of the ratio spectrum R(v) (averaged between epochs) in 
every velocity bin from -30,000 to km s _1 for all the 
sources in our sample. R a ii roughly represents the dis- 
tribution of absorption depths in all velocity bins of the 
"average BAL" in our sample. 

We then compare R a ii to the distribution of the ra- 
tio spectra for the more-absorbed epoch and also for the 
less-absorbed epoch. We define (il<), as the ratio spec- 
trum, R(y), averaged over the bins in a varying region in 
the epoch for which absorption is stronger (and therefore 
R(v) is smaller). We take the more-absorbed epoch for 
each individual region, and for some QSOs with multiple 
regions (e.g., 1208+1535), this means we take the average 
ratio spectrum in the first epoch for some regions and in 
the second epoch for other regions. Similarly, we define 
(i?>), obtained from the epoch with weaker absorption 
(and therefore a larger average R(v)) for each region. 

In Figure [7) we plot the distributions of (R<) (top 
panel), (i?>) (middle panel), and R a u (bottom panel, 
with arbitrary normalization). The distribution of ab- 
sorption depths in the varying regions in the epoch of 
weaker absorption ((i?>)) resembles that of the average 
BAL absorption (R a ii), while the distribution from the 
more strongly- absorbed epoch ((il<)) differs consider- 
ably from Rail- This may simply be a consequence of 
the fact that most bins in the 1400-1550 A region are 
not strongly absorbed for the average BAL, and the dis- 
tribution il>, taken from the less-absorbed epoch, more 
accurately reflects this. 

The average velocity of varying regions with (i?>) > 
0.9 is -19,000 km s _1 , compared to -13,000 km s _1 for 
those with (i?>) < 0.9. A Kolmogorov-Smirnov (KS) 
test suggests (at 99.3% confidence) that the velocities for 
these two sets of varying regions are distributed differ- 
ently. It is not surprising that varying regions which are 
weakly absorbed in one epoch tend to appear at higher 
velocities, as BAL abs orption is, on aver age, weaker at 
higher velocities (e.g.. iKorista et aLlll993l ). 

A total of 20 of the 30 varying regions have (R<) > 0.5, 
indicating a preference for measurable variation to occur 
in regions with weaker absorption. Half of the varying 
regions have (i?>) > 0.9, corresponding to weak or even 
no absorption in one epoch. In the (R a ii) distribution, 
88% of bins have (R) > 0.5, and 57% have (R) > 0.9. 
The tendency for variation to occur in weakly-absorbed 
regions is therefore a reflection of the fact that most bins 
in our sample spectra are weakly absorbed, combined 
with the tendency for variation to occur at a wide range 
of outflow velocities. We find no evidence for a correla- 
tion between the strength of variation in a varying region, 
| (Ail) |, and the average depth of that region, (R). We 



therefore find no evidence that the probability of varia- 
tion in a velocity region is dependent on the absorption 
strength in that region. 

3.1.4. Fraction of Absorption Width That Varies 

We define Jbal as the fraction of the 1 A wide wave- 
length bins from 1400-1550 A which are at least 20% 
absorbed (i?(A) < 0.8) in one or both epochs. We de- 
termine the fraction of varying wavelength bins, f vo , by 
counting the number of 1 A bins from 1400-1500 A which 
vary by at least |Ai?(A)| =0.2 between epochs. (Here we 
require a greater threshold of change |Ai?(A)| than in 
N3.1.1I because we are considering single bins.) Because 
all spectra are binned onto the same 1 A grid which has 
149 bins, the fractions Jbal and f vo are both propor- 
tional to the number of deep or varying 1 A bins, with 
the same constant of proportionality, 149. We do not cap 
the ratio spectra at 1 when determining /„&, as this could 
bias our results. We have not accounted for absorption 
from intervening systems, so QSOs with strong interven- 
ing system lines in the C IV region may have a slightly 
higher Jbal than measured from the BAL alone. 

A Spearman rank correlation test finds that /„{, and 
Jbal are correlated at 99.9% confidence. That is, BALs 
with broader regions of deep absorption (-R(A) < 0.8) 
tend to have a larger number of variable wavelength bins 
from rest to 30,000 km s _1 . A linear fit gives: 

fvb = -0.01 ± 0.03 + (0.39 ± 0.07) f bal- (1) 

The errors are estimated from the scatter in the data 
(rather than from unknown error on indi vidual data 
points ), assuming xt = 1; see pp. 666-669 of lPress et al.1 
(2002). Figure [8] shows the plot of f vo against Jbal for 
our sources together with the fit from Equation [TJ 

Jbal and Sa« (defined in ^3.1.1|) are correlated at the 
98% confidence level, supporting the finding that broader 
(deep) BALs have more varying absorption bins. If we 
confine the search for varying bins to only include those 
which are deeply absorbed (i?(A) < 0.8), we find a simi- 
larly strong correlation (at 99.8% confidence) as varying 
bins primarily occur in the strongly-absorbed regions. 

3.2. EW Variation 

In §3.1[ we determined ranges of velocity bins over 
which the BAL was observed to vary beyond a certain 
threshhold. The strength of the variation in a particu- 
lar velocity bin was not important, so long as it exceeded 
the threshold. By contrast, tests for variation in EW also 
account for changes in absorption depth. We calculate 
EW errors formally, estimating a continuum error of 5%. 

3.2.1. Cumulative EW Variation 

In Figure [9l we plot the total change in equivalent 
width between epochs, AEW, against the average EW, 
(EW), for the entire C IV region of each BAL QSO in 
our sample. For comparison, we have also plotted the 
results from the studies of B93 and L07, which apply to 
shorter time scales. In the four cases where L07 mea- 
sured EWs from two BALs in a single QSO separately, 
we have combined the EWs to get an approximate EW 
for the entire C IV absorption region. In cases where B93 
observed a source more than twice, we used data from 
the two epochs with maximal time separation. 
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We do not find significant correlations between (EW) 
and AEW in the combined sample. There is a signifi- 
cant (>99.99% confidence) correlation in the combined 
sample between (EW) and the magnitude of fractional 
change in EW, \AEW/ (EW)\, so that the magnitude of 
fractional variability is greater for weaker BALs. How- 
ever, this correlation is primarily driven by the weakest 
BALs. It disappears (confidence < 90%) if we exclude 
the 20% of BALs (mostly taken from L07) in the com- 
bined sample which have EWs >-1000 km s _1 . For this 
reason, we do not draw strong conclusions from the cor- 
relation. 

We plot AEW against the logarithm of the change in 
flux density at 1400 and at 2500 A in Figure Q3H We 
find no significant correlation between the change in flux 
density, AF V , and AEW, or between \AF V \ and \AEW\, 
as might be expected in a photoionization-induced vari- 
ation scenario. Ionization fractions for a given ion are 
roughly symmetric about a peak at some ionizing flux 
(for a given shape of the ionizing continuum). While the 
ionization fraction for a given ion can remain unchanged 
if the ionizing flux changes by just the right amount, this 
is a contrived scenario, and is unlikely to occur in a large 
number of cases. We discuss physical causes of EW vari- 
ation further in £14.31 and £14.41 

The bottom panel of Figure fTOl shows that the changes 
in F„ at 1400 and 2500 A are correlated. The photoion- 
ization of C IV is driven by photons at higher energies 
(>50 eV) than we have observed. If the high-energy con- 
tinuum drives photoionization variation, it would have 
to vary independently of the observed continuum (at 
A > 1400 A). The nature of QSO continuum variation 
from the UV to X-rays is not currently well-constrained; 
however, this phenomenon has received some attention 
for lower-luminosity Seyfert galaxies. In the case of 
the Seyfert 1 galaxy NGC 3516, the X-ray and opti- 
cal continua have been observed to vary independently 
(|Maoz et al.l l2002f ) . Studies of other Seyfert 1 galaxies 
have found time lags of tens of minutes to days between 
the X-ray and optical con tinua (e.g., IShemmer et al.l 
120031 iDesroches et all 1200^ 1. 

3.2.2. EW Variation at Different Outflow Velocities 

In order to test for any velocity-dependence in BAL 
variation, we calculate the EW for each BAL in 
5000 km s^ 1 wide bins from -30,000 to km s" 1 . We 
find no evidence of correlation between At sys and either 
AEW or \AEW/(EW) | in any of the velocity bins. How- 
ever, the range of At sys in our sample covers only a fac- 
tor of ~2, so we are not sensitive to velocity-dependent 
variation on a wide range of time scales, nor are we sen- 
sitive to faster variation which does not generate trends 
on multi-year time scales. 

BAL absorption is likely saturated with a depth 
primarily determ ined by the covering factor (e.g., 
lArav et al]|1999al) . BAL absorption is, on average, shal- 
lower at high velocities. If shallower absorption re- 
gions are less saturated, high-velocity absorption could 
be more responsive to changes in the BAL absorber. 
However, we find no indication (either visually or for- 
mally with a KS test) that the distribution of AEW 
changes across velocity ranges. 

Similarly, if absorption at higher velocities were less 
saturated, high-velocity variation may track the contin- 



uum more closely. We test for correlations between AF U 
at 2500 A and AEW in 5000 km s" 1 velocity bins. The 
strongest correlation is between |A£W| and AF„ in the 
-10,000 to -15,000 km s" 1 velocity range. At 99% confi- 
dence, the putative correlation is not highly significant, 
given that we have tested 12 cases for correlations (both 
AEW and |A_EIV| against F v in 6 velocity bins). We 
therefore do not find any strong evidence for velocity- 
dependent correlations with the continuum. 

3.2.3. Evolution of \ AEW \ 

In order to compare the L07 sample to the BALs in 
our study, we construct a matched sample of L07 objects 
which have EWs greater in magnitude than that of the 
weakest BAL in our sample. We call this set of 18 objects 
the "strong L07 BALs." Figure QT] shows the change in 
EW, AEW, plotted against the time between epochs. 
The data from B93 and the strong L07 BALs are in- 
cluded. The range (or "envelope" ) of AEW is clearly in- 
creasing with time. The B93 BALs were slightly stronger 
on average than our BALs, but had smaller |Ai?W| over 
shorter time scales. For the strong L07 BALs, the mean 
|AJ5W| of 390 km s~ x is much smaller than that of our 
sample (1300 km s _1 ). 

In Figure [TTJ we plot the standard deviation of a slid- 
ing window of 15 time-ordered objects to illustrate the 
increasing spread of AEW. The mean value of |A_EW 
in our sample is 1300 km s" 1 over a mean time of 
4.3 yr, while the standard deviation is «1650 km s _1 
at At sys w 4 yr. 

A strong BAL with (EW) = -6000 km s" 1 which 
varies at the rate of 400 km s -1 yr~ x would have a char- 
acteristic lifetime of \(EW) / (AEW/ At sys )\ « 15 yr. If 
the envelope surrounding AEW keeps growing on longer 
time scales, some strong BALs could completely disap- 
pear over a few decades. 

We might also expect that strong BALs would appear 
in QSOs that had previously shown no BALs. (We dis- 
cuss constraints on BAL transience further in ^4.20 Fu- 
ture studies with a large number of BAL observations 
over a wide range of time scales will allow careful study 
of the time-evolving distribution of AEW, which is ap- 
parently revealing important physics of the QSO envi- 
ronment. 

3.2.4. Evolution of \AEW / (EW)\ 

BALs have been observed to vary on short (multi- 
month) time scales §1.3|) . and L07 found that the 
largest fractional variability \AEW/(EW)\ was seen in 
sources with the longest (« 4 month) rest-frame time be- 
tween epochs. If BAL QSOs evolve monotonically over 
several years, we would expect to see larger variation on 
the longer time scales in our study. 

The mean value of \AEW/(EW) \ is 0.33 ± 0.04 for all 
L07 objects and 0.32 ±0.04 for our sample. However, the 
mean of the L07 sample is dominated by the three outlier 
sources with \ AEW/(EW)\ ~ 1.5, and all three outliers 
are from sources with weaker BALs than in our sam- 
ple. For the strong L07 BALs alone ( £J3.2.3|) . the mean 
is 0.12 ± 0.02. For the B93 sources, the mean fractional 
change in EW is only about 0.12±0.06, similar to that of 
the L07 sources with stronger BALs. A KS test indicates 
that \AEW/(EW) \ is distributed differently between the 
B93 sources and our sources at >99.9% confidence. The 
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distribution of fractional change in EW also varies be- 
tween the strong L07 BALs and our sample at 99.3% 
confidence. 

We have plotted our values of \AEW/(EW)\ against 
time in Figure [12] along with a matched sample of values 
from B93 and L07. We have only included sources with 
| (EW) | greater than the minimum of our sample in the 
plot. 

In summary, we find that (apart from a few outliers) 
the variation in fractional EW is greatest over the longest 
times. For three outliers with relatively small (EW) in 
the L07 sample, \AEW/(EW)\ has varied greatly over 
short (2-4 month) times. 

3.3. Acceleration at Detachment Velocity 

The relatively long time scales in our study al- 
low a sensitive search for acceleration of BAL com- 
ponents. Previous studies have found accelerations 
of sa 0.03 cm s~ 2 for Q 1303+308 (IVilkoviskii fc Irwin 
2001h a nd sa0.08 c m s~ 2 for Mrk 231 (jRupke et al 
20021 ) . lHall et alJ (|2007h have recently presented ev- 
idence for acceleration in a C IV BAL trough of 
SDSS J024221.87+004912.6 at a rate of w0.15 cm s" 2 
over 1.4 yr. The detachment velocity of this BAL trough 
increased in magnitude between epochs by «70 km s . 
The wavelength calib ration of the LBQ S spectra is ac- 
curate to about 4 A (|Foltz et al.lll987al) . corresponding 
to about 260 km s _1 at 1549 A for a z = 2 QSO. If the 
BALs in our sample accelerate over 6 yr at the same rate 
as seen in SDSS J024221. 87+004912.6, we should detect 
noticable trough changes. 

The variation in SDSS J024221.87 + 004912.6 was most 
evident in the onset region of the deep C IV BAL. The 
objects in our sample with deep BAL troughs generally 
have steep onset regions (at sa-5000 km s _1 ) which are 
ideal for tests of acceleration (Figure [3]) . After convolv- 
ing the SDSS spectrum with a Gaussian to approximate 
the LBQS resolution, we do not find strong candidates 
for acceleration in the BAL onset region. This does not 
necessarily mean that material is not accelerating along 
the outflow; the onset region may be continually replen- 
ished with material at a constant velocity. We discuss 
the physical implications of this measurement in 

The ratio spectrum, R(X), drops below a threshold of 
i?(A) =0.3 at the same wavelength (within 1 A) in both 
epochs for the 7 sources listed in Table [5j The upper 
limits on acceleration (Av / At sys ) for the seven sources 
with sharp BAL onsets range from 0.12 to 0.17 cm s~ 2 
and are given in the table. The acceleration limits were 
calculated using the different rest-frame time elapsed be- 
tween epochs for each source, assuming an upper limit of 
1 A for movement of the absorption region. In 3 of the 
7 sources (0018+0047, 0051-0019, and 1331-0108), nar- 
row absorption lines from apparent intervening systems 
are well-matched between epochs and further support the 
use of 1 A as a constraint on the relative wavelength cal- 
ibration. 5 The velocities of the intervening systems are 
given in Table [2j 

5 Shifting the narrow absorption regions of the SDSS ratio spec- 
trum in wavelength space and fitting against the LBQS ratio spec- 
trum indicate that the wavelength calibration for these sources 
may be good to 0.5 A or less. However, unmodeled emission and 
systematics prevent us from using the x 2 statistic to place strong 
confidence limits on the relative wavelength calibration. 



An eighth source, 1208+1535, also has a steep onset 
region, but the BAL variation is so complex that it is 
difficult to determine whether acceleration is a factor. 
Shifting the later (SDSS) epoch spectrum to lower veloc- 
ities by about 2500 km s^ 1 would match some features in 
the LBQS-era BAL trough, but other features would not 
match well. This QSO is an excellent target for future 
monitoring. 

The variation in most absorbed regions of the sources 
in our sample (excepting the onset region of strong 
BALs) is too complex to test for acceleration. Visual 
inspection of weaker absorption regions shows some vari- 
ation between epochs, including narrowing of the absorp- 
tion region (e.g., in 0010-0012 at ^-13,000 km s" 1 ) or 
variation on only one side of the trough (e.g., 1234+0122 
at ~-13,000 km s _1 ). Only the lowest-velocity compo- 
nent of 0109-0128 (at «-3000 km s" 1 ) shows visual ev- 
idence of acceleration on both sides of the absorption 
trough. Unfortunately, in this case there are no narrow 
lines from intervening or other intrinsic systems strongly 
present in both epochs which can be used to test the 
relative wavelength calibration. The doublet structure 
of the feature (seen in the full resolution of the SDSS 
spectra) indicates that the feature, if intrinsic, could be 
attributed to C IV flowing outward at -3400 km s _1 . We 
do not see any clear cases of significant absorption from 
lines of other ions at this outflow velocity in either the 
LBQS or SDSS spectrum. Future spectroscopy of this 
object will be able to determine whether or not this is 
an accelerating, intrinsic absorption component. If this 
feature is actually accelerating, and continues to acceler- 
ate at the same rate, it will join the main BAL trough in 
about 70 yr (in the observed frame). 

4. DISCUSSION 

In the following sections, we consider the physical im- 
plications of our study. First, we present a brief sum- 
mary of the main quantitative results of our analysis 
in ij4.ll We then briefly discuss constraints that our 
study places on BAL transience (jj4|2]). We also consider 
the implications our study has for the physical processes 
of photoionization-induced variation fi )4.3p and geomet- 
ric variation ( §4.4)1 . Finally, we discuss the implications 
of our acceleration constraints on the BAL onset region 

4.1. Summary of Variation Metrics 

Although C IV BAL absorption can extend up to (and 
even beyond) -30,000 km s _1 and BAL troughs can be 
tens of thousands of km s _1 wide, variation tends to oc- 
cur in (multiple) narrower bands of width <5000 km s _1 , 
with smaller widths (<2000 km s^ 1 ) most common. The 
varying regions of all sources combined are distributed 
across a wide range of velocities. We find no evidence 
that the velocity width of a varying region is correlated 
with the outflow velocity of that region. 

The degree of variation, measured in terms of the av- 
erage change in the ratio spectrum across the varying 
regions, |(Ai?)|, peaks at \(AR)\ ~ 0.2. The algorithm 
used to determine varying regions generally requires that 
| (Ail) | > 0.1 (set by the data quality), so we are not able 
to identify weak, broad variation. The degree of varia- 
tion | (Ail) | is not strongly correlated with the absorbed 
depth, velocity, or velocity width of outflowing regions, 
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nor with the elapsed time between observing epochs. 

About half of the varying regions are associated with 
components which are seen strongly in only one of the 
two epochs. These "transient" regions tend to appear 
at higher outflow velocities (-19,000 km s _1 ) than the 
non-transient regions (-13,000 km s" 1 ). This is likely a 
consequence of the fact that BAL absorption tends to 
be shallower at higher velocities, while variation occurs 
across the entire BAL trough. 

The number of bins which vary in the C IV absorption 
region is correlated with the number of deeply-absorbed 
bins. The average value of \AEW/(EW)\ is » 0.3 for 
our sources, somewhat larger than for a matched sam- 
ple of L07 and B93 sources observed over shorter time 
scales. We do not find evidence that EWs vary differently 
at different velocities, even though the ave rage BAL ab- 
sorpt ion varies strongly with velocity (e.g. jKorista et all 
1993). As Fi gure [Tl] clearly shows, the distribution of 
AEW broadens with time. 

We do not find any strong evidence of BAL absorber 
acceleration in our sample up to wavelength calibration 
limits, although the variation was so complex in several 
cases that it could have disguised acceleration. The steep 
absorption dropoff at the detachment velocity did not 
appear to change velocity between epochs. We found 
upper limits to acceleration of the onset region of 0.12 
to 0.17 cm s~ 2 . The absorbing material in the onset 
region is being replenished in such a way as to generally 
preserve the trough shape, is accelerating intermittently, 
or is generally accelerating less rapidly than in the case of 
SDSS J024221. 87+004912.6. We did identify ambiguous 
cases for acceleration which we could not constrain, but 
will test with future observations. 

4.2. BAL Transience 

The degree of variation observed in our study raises 
the possibility that some strong BALs could be transient 
on observable time scales f ^3.2.3[) . We briefly present 
constraints on BAL lifetimes in this section, although 
the limited nature of our sample, both in numbers and in 
time scale, prevents us from drawing strong constraints. 
Future studies will obtain improved results from large 
samples of BAL QSOs in surveys like SDSS. 

Suppose BAL outflows are a common, but transient, 
phenomenon in QSOs. If a given QSO with lifetime tQso 
hosts one BAL covering the entire QSO sky with a life- 
time tBAL <C tQso, and if we take two observations of 
the QSO separated a time t i, s <C tBAL apart, then we 
have a chance rj t \, s /tQso of catching the BAL as it 
appears, assuming the BAL forms quickly. If a QSO 
hosts Nbal such BALs over its lifetime (but at most 
one BAL at a given time), we can detect a BAL ap- 
pearing during a fraction j3 = NsALtobs/tQso of the 
QSO lifetime. In a simplified evolutionary scenario in 
which each QSO hosts BALs for 20% of its lifetime (based 
on the frequency of B ALs in QSOs; iBecker et alj|2000t 
iHewett fc Foh9 120031 ). we have N B ALt B AL = 0-2t QS o- 
In that case, (3 = 0.2t o b s /tBAL- We observed 115 QSOs 
with no significant BALs in either epoch. If BALs are 
equally probable to appear at any time, then (3 is also 
the probability that we observe a BAL forming in a 
given QSO, and a binomial distribution indicates there 
is a > 90% chance we would have seen one BAL appear 
out of 115 chances if /? > 0.02. In this scenario, the 



BAL lifetime is tsAL > 43 yr and the QSO lifetime is 
tQSO > 210N BAL yr. 

If a strong BAL has a lifetime tBAL, we would ex- 
pect BALs observed in the first epoch to disappear in 
a fraction t b s /tBAL of the observations. Given 9 QSOs 
determined to have strong BALs (Bio > 100) during a 
first observation epoch, if we make a second observation 
tobs ~4.3 yr later and the BAL lifetime is <18 yr, we 
have a >90% chance of observing at least one case in 
which a strong BAL disappears (using a binomial distri- 
bution as before). We therefore expect tsAL > 18 yr. 
BAL lifetimes may, of course, be at least several orders 
of magnitude longer. 

4.3. Does Photoionization Variation Drive BAL 
Variation ? 

Although BAL profiles are likely primarily determined 
by the velocity-dependent geometric coverin g factors of 
saturated absorbers (e.g.. lArav et al.lll999af) . we briefly 
consider a hypothesis where the trough depth is at least 
partly determined by the (unsaturated) absorber column 
density. In this case, the BAL trough could vary in re- 
sponse to the ionizing continuum. But, as in previous 
studies, we have not found any significant correlations 
between (the magnitude of) absorber variation and con- 
tinuum variation. We have tested for correlations with 
the continuum at both 1400 and 2500 A. Of course, the 
C IV ionization state will be strongly influenced by the 
continuum at shorter wavelengths where optical studies 
do not have spectral coverage. While the continuum vari- 
ations at 1400 and 2500 A are highly correlated, the far 
UV and X-ray continuua may vary independently of the 
optical continuum. As discussed in £13.2. 1[ this has been 
observed to happen in at least one Seyfert 1 galaxy. 

We find no significant evidence that absorption vari- 
ation is dependent on absorption depth (i )3.1.3p . as we 
might expect if weakly-absorbed troughs tended to be 
less saturated (and thus more responsive to changes in 
ionization levels). In fact, we observe specific cases where 
deeply-absorbed regions vary, or contiguous regions at 
similar absorption depths vary independently (e.g., 0051- 
0019). 

The broadening distribution of |A£W| with time 
( t|3.2.3[) and the lack of variation in some parts of BALs 
is not likely to be due simply to long photoionization 
time scales. The gas density (estimated from the C IV 
recombination time scale) would need to be <7000 cm" 3 
to have response times >4.2 yr. For typical equilibrium 
photoionization models, this would require the BAL ab- 
sorber to be at least hundreds of parsecs from the ionizing 
source. In comparison, the broad emission line regions 
of QSOs are believed to be less than a few tenths of a 
parse c from the central source (e.g., iKaspi et all [2000. 
120071 and references therein). 

4.4. Covering Factor Variation 

The observed absorption in BAL troughs is strongly in- 
fluenc ed by the absorber covering factor fe.g. JArav et al.l 
fl999a|) . In this scenario, the absorption may be satu- 
rated, but if the absorber does not cover the continuum 
source completely, the absorption is not black. In a gen- 
eral sense, covering factor variation also includes scenar- 
ios w here the continuu m is scattered around the absorber 
fe.g.. lOgle etaUH999h . 
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Several factors in our study support a scenario in which 
the absorber geometry varies. BALs in our sample vary 
in discrete velocity segments up to about 5000 km s _1 
wide (but typically < 2000 km s _1 wide) which are nar- 
row compared to the extent of a strong BAL. The change 
in absorption depth is commonly |(Ai?)| «0.2 in a vary- 
ing region. This leads us to speculate that the BAL ab- 
sorber contains clumps of material which extend across 
<2000 km s _1 in radial velocity and cover «20% of the 
continuum emitter (assuming the clumps are optically 
thick). 

Seen individually these clumps may appear similar to 
narrower absorption features such as mini-BALs or even 
NALs. It has been previously suggested that BAL ab- 
sorption may be an extension of th e NAL phenomenon 
(e.g., §4.2.1 of lArav et all fl999bh . and narrower fea- 
tures, includi ng C IV doublets, have been observed 
to vary (e.g.. lHamann et all 119951: iGangulv et alj 120011 : 
INaravanan et al.l l2004t I Wise et alj 120041) . Modeling is 
needed to determine if and how it is possible to gen- 
erate clumps of material which are dynamically linked 
over (up to) 2000 km s _1 in the BAL outflow. Previ- 
ous observations of absorption variation have determined 
that such structures apparently can exist. For example, 
lHamann et alj (|1997| ) observed intrinsic N V, Si IV, and 
C IV absorbers varying in unison over <4 months in the 
QSO Q2343 + 125. Each of the C IV doublet compo- 
nents was «400 km s _1 wide (FWHM), and the overall 
varying absorption feature was wlOOO km s" 1 wide. The 
absorbers were found to cover <20% of the continuum 
emission, and there was some indication that the varia- 
tion was driven by c hanges in covering factor. Recently, 
iMisawa et al.l (|2007| ) observed a mini-BAL with signifi- 
cant substructure that varied in concert in the QSO HS 
1603+3820. 

Over 3-6 years, about one third of a BAL varies. Apart 
from a few outliers, AEW/ (EW) increases with time; 
it is «0.3 on average in our sample. These results are 
also similar to variation properties of narrower absorp- 
tion features, as at least 20-25% of associated absorp- 
tion lines (AALs; narrow lines within 5000 km s -1 of the 
emitter rest frame) va r y on multi-year t ime scales (e.g., 
INaravanan et al.ll2004t (Wise et al.ll2004f ) . 

Several specific cases of BAL variation arc qualita- 
tively similar to what we would expect from covering 
factor variation. We see cases (most notably in 0109- 
0128) where the change in absorption strength, Ai?(A), 
is nearly constant across a wide velocity range. This 
could occur in the limiting case where the covering fac- 
tor is independent of the outflow velocity. Wc also sec 
cases (most notably in 0051-0019) where only part of a 
wide region (all at a common absorption strength R(X}) 
varies. This could occur if absorbing clumps at some ve- 
locities moved out of the line of sight, while clumps of 
the same size (i.e., covering factor) remained in view at 
other velocities. Of course, it is possible to manufacture 
covering factor scenarios for arbitrary patterns of BAL 
variation. Larger samples of BAL variation on a range of 
time scales would improve our understanding of the gen- 
eral characteristics of BAL variation and would address 
the plausibility of such scenarios. 

Studies with sufficient spectral resolution to resolve 
doublets in both epochs may also be able to constrain 
covering factors quantitatively, at least in cases where 



the absorber has doublet s ubstructures. The mini- 
BAL recently monitored by IMisawa et al.l ((2007) con- 
tained several narrow sub-components, yet the entire 
mini-BAL varied in concert. In this case, the varia- 
tion was attributable particularly to changes in cover- 
ing factor. In another recent study, I Hall et aTl (|2007l ) 
concluded from an analysis of the doublet structure of 
Si IV absorption components that the BAL absorbers in 
SDSS J024221. 87+004912.6 were composed of multiple 
clouds or filaments with scale sizes ~10 15 cm. 

The fact that the envelope for AEW increases with 
time (Figure ITTj) indicates that BAL outflows are chang- 
ing along the line of sight on observable time scales. If the 
variation is caused by geometric factors such as clumps 
moving across the line of sight, the pattern of varia- 
tion is revealing information about the distribution of 
clumped material in the outflow. The clumps are appar- 
ently common in BALs at velocities ranging from -5000 
to -25,000 km s _1 . Because we see variation in even the 
deepest absorption regions, the clumps may not be so nu- 
merous at any velocity as to obscure the source several 
times over (at least in some cases of deep absorption). 

4.5. Acceleration 

We searched for acceleration at the onset velocity of 
strong absorption features previously observed in 7 BALs 
but found no evidence for such acceleration in our sam- 
ple. Over 3-6 yr, acceleration due to radiation pres- 
sure on a C IV absorber would certainly be noticable in 
the sharp absorption onset regions in our sample, unless 
the absorbers are far enough (parsecs or more) from the 
emitter to geometrically dilute the radiation. It has been 
suggested that the onset region is being replenished (at 
constant velocity) from the accretion disk, so that the ab- 
sorption onset velocity does not cha nge even though ab- 
sorber material is accelerating (e.g.. iMurrav et al .1119951 : 
iProga et a l. 2000). In this case, the onset region may 
represent the velocity at which the wind leaving the disk 
is bent into our line of sight by radiation pressure. 

Similarly, the tight upper limit on acceleration of the 
C IV mini-BAL in th e Seyfert 1 galaxy NGC 4151 led 
IWevmann et al.l (|1997l ) to conclude that the absorber was 
either very distant from the central source, was slowed 
by drag in a surrounding medium, or had a complex ge- 
ometric structure. A small acceleration upper limit of 
0.1 cm s~ 2 was also found for li nes from several i ons in 
the Seyfert 1 galaxy NGC 4051 rtKaspi et al.ll200l . 

Over the course of our observations, the accretion disk 
has rotated significantly. For a 10 8 M Q black hole, the 
inner edge of the BAL wind is esti mated to be at a radius 
of ~10 16 cm (jMurrav et al.lll995h . Material in a Keple- 
rian orbit at 10 16 cm (700 gravitational radii) from such 
a black hole would complete more than two revolutions 
in 4 yr. A wind leaving the disk would presumably share 
the disk rotation. Given the observed (BAL) structure 
in the accelerated outflow, the possibility of structure in 
the accretion disk itself, and variation in the illuminat- 
ing source, it is remarkable that the onset regions in our 
sample are so stable. At a radius of 10 17 cm, the amount 
of rotation would be only about 8% of one revolution. 

5. SUMMARY AND FUTURE WORK 

We have searched for C IV BAL variability in QSOs 
which were observed in both the LBQS and SDSS (or by 
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the HET) . These surveys were conducted approximately 
two decades apart (in the observed frame) , enabling com- 
parisons between QSOs at 1.7 < z < 2.8 over rest-frame 
time scales of 3-6 yr. Because our sources were (with one 
exception) drawn from two large-scale QSO surveys, we 
expect that they constitute a reasonably representative 
sample of optically-selected BAL QSOs in the allowed 
rcdshift range. We summarize our most important find- 
ings here. 

1. BALs tend to vary on multi-year time scales in ve- 
locity regions which are a few thousand km s _1 
wide, much narrower than a strong BAL. The vary- 
ing regions occur at a wide range of outflow veloc- 
ities and absorption depths. The typical change in 
absorption depth is <25% of the continuum. 

2. The number of spectral bins which vary over 3-6 yr 
is correlated with the number of deeply-absorbed 
bins, and the magnitude of the fractional change 
in EW, \AEW/(EW)\, is about 0.3 on these time 
scales. 

3. The range of AEW for BAL QSOs increases with 
time, up to 3-6 yr. Because we did not see strong 
BALs appear or disappear, we constrain typical 
BAL lifetimes to be at least 18 yr. However, if the 
envelope of AEW continues to increase on longer 
time scales, some strong BALs could have lifetimes 
as short as a few decades. 

4. We find no evidence that the variation is dominated 
by photoionization on multi-year time scales. The 
variation does not correlate with changes in the 
observed continuum, although we note that C IV 
ionization levels are driven by photons at higher en- 
ergies than we observe. Even if the ionizing contin- 
uum varies independently of the optical/UV con- 
tinuum, the patterns of variation we observe would 
require complicated outflow structures. 

5. Several aspects of the variation we observe are 
at least qualitatively consistent with a scenario in 
which the covering factor is changing on multi-year 
time scales for clumps of material in the BAL out- 
flow. 

6. We find no evidence for acceleration in the steep 
onset region of BAL troughs, despite the expected 
radiation pressure on the absorber. In scenarios 
where the onset region is replenished by material 
flowing off the accretion disk, the disk would have 
rotated significantly over 4 yr unless the outflow is 
at a radius > 10 17 cm (for a 10 s M Q black hole). 
The flow would have to be remarkably azimuthally 
smooth not to show evidence of this rotation. 

Previous multi-wavelength studies have searched for 
correlations wi t h BA L properties. For example, 
iGallagher et al.l (|2006l ) found a correlation between the 
maximum BAL outflow velocity and the X-ray weakness 
for a sample of 35 LBQS BAL QSOs, but did not find evi- 
dence for other correlations. In multi- wavelength studies, 
the different wavebands are typically observed some time 
apart. Our study shows that BALs can evolve signifi- 
cantly over typical time scales between observations, in- 
troducing additional scatter into relations between BAL 
and, e.g., X-ray properties. 



11 

Together with the previous work of B93 and L07, we 
have shown that studies of BAL variability on human 
time scales yield interesting and important results. The 
future of such studies is promising, as the number of 
known BAL QSOs has incre ased dramatically w ith the 
SDSS. Based on the results of lTrump et al.l (|2006f ). we es- 
timate that the SDSS DR5 QSO catalog contains «3000 
sources at z > 1.7 with measurable, traditional C IV 
BALs (having BI > 0) and «7500 with C IV absorp- 
tion features at least 1000 km s _1 wide. Observations of 
these sources with large-scale, spectroscopic surveys will 
an increase in the size of BAL variability surveys by or- 
ders of magnitude, and will enable more sensitive studies 
of evolution on a range of time scales. 

QSOs must be significantly redshifted in order to move 
the C IV absorption region into the optical spectrum. At 
z = 2, covering a time range only three times longer than 
the 3-6 yr span in our study would take nearly a typical 
professional lifetime (27-54 yr). It will therefore be dif- 
ficult to increase the time scale for C IV BAL variation 
studies greatly unless UV spectra of lower-redshift QSOs 
are obtained. There is also much to be gained by improv- 
ing the coverage on time scales <10 yr. The distribution 
of AEW and the relation between BALs and narrower 
absorption features will be interesting to map out for a 
range of time scales and intrinsic QSO properties. In- 
creased spectral resolution will also assist in determining 
absorber substructure and constraining covering factors. 
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APPENDIX 
NOTES ON INDIVIDUAL SOURCES 

In this appendix, we briefly discuss the spectra of the C IV regions for the 13 BAL QSOs in our study. In order to 
compare C IV absorption spectra from two epochs directly, we constructed the ratio spectrum, i?(A), for each epoch 
by dividing out the continuum and broad line emission as described in §2.11 We convolved the SDSS spectrum with a 
Gaussian of the appropriate width to approximate the LBQS spectral resolution. The ratio spectra for the BAL QSOs 
in our sample are shown in Figure [3l 

We have visually searched the SDSS spectra for absorption from narrow doublets up to w 100, 000 km s _1 blueward 
of the broad Mg II A2799 emission line. These doublets are likely due to Mg II ions in interveni ng absorbers. After 
determining the velocity offset of these (likely) intervening systems, we search a theoretical line list (|Verner et al.lll996[) 
for lines which, at that velocity, may also appear in the C IV absorption region. We do not know the ionization state 
or elemental abundances of the intervening systems, so we cannot predict the strengths of these lines. However, their 
velocity widths and equivalent widths are small compared to intrinsic BAL absorption features, and most of these lines 
do not reside in variable regions anyway (Figure [3]), so they will not significantly influence our results. The locations 
of (potential) interesting lines with absorption oscillator strengths > 0.1 arc marked in Figure [3l The velocities of 
these intervening absorbers are given in Table [5J 

0010-0012 

A relatively narrow BAL has decreased in width over time. BI = in the SDSS observation. 

0018+0047 

This source has a strong BAL at a relatively low outflow velocity (-2000 to -6000 km s _1 ). The absorption has 
perhaps weakened between epochs. Intervening absorption may be weakly present in the C IV absorption trough. 

0021-0213 

In this case, the second epoch is an HET LRS spectrum, rather than an SDSS spectrum. Weak variation is seen in 
narrow components, while a stronger trough has appeared at -16,000 km s _1 . 

0051-0019 

Four strong, broad C IV absorption troughs are visible at velocities ranging from -5000 to -22,000 km s _1 . The 
weaker C IV absorption components show the greatest variation. The low-velocity region of the highest-velocity 
absorption component (at -18,000 km s _1 ) has remained nearly constant while the remainder of the high-velocity 
component has disappeared. The broad absorption at w-18,000 km s _1 may, in principle, be contaminated by a 
narrow line from Al I in an intervening system, but the full resolution of the SDSS spectrum shows that the absorption 
at that velocity is too broad (~2000 km s _1 wide) for the contamination to have any significant effect. 

0055+0025 

This source is technically classified as a BAL QSO based on weak, broad C IV absorption in the -11,000 to 
-14,000 km s -1 range. The absorption feature is weak, and disappeared between epochs. 

0109-0128 

This source has varied dramatically between epochs. Two broad, high- velocity components have appeared at -23,000 
and -17,000 km s _1 . Weak absorption features may be present at lower velocities in the LBQS epoch. The broad, 
low-velocity BAL between -12,000 and -5,000 km s _1 has strengthened considerably, almost to black at the deepest 
point. Between -10,000 and -6,000 km s _1 , the difference between ratio spectra is nearly constant, consistent with 
a scenario in which the covering factor of the source varies for a saturated absorber. The low-velocity absorption 
component at -3,000 km s _1 shows evidence of acceleration between epochs, but the putative acceleration (of magnitude 
w200 km s _1 ) is smaller than the accuracy of the wavelength calibration of the LBQS spectrum. 

1208+1535 

Dramatic variation is evident in this case. The deepest points in the C IV absorption and the local maximum 
between them have shifted several thousand km s _1 to higher velocities. The onset velocity (about -7000 km s _1 ) has 
not changed. Broad, shallow absorption has increased from -17,000 to -23,000 km s _1 . 

1213+0922 

Technically classified as a BAL QSO based on our BI measurements, 1213+0922 was not classified as an HF03 BAL 
QSO. The absorption is shallow and broad. We note that similar ratio spectrum shapes are seen in several SDSS QSOs 
(not observed in the LBQS and therefore not in our sample). We cannot rule out the possibility that the feature is 
due to structure in the continuum emission, rather than broad absorption. (See also £)2.2I ) 



BAL Variation on Multi-Year Time Scales 



13 



1234+0122 

This source was classified as an HF03 BAL QSO, but our continuum fit technically gives BI = in both epochs. 
Nonetheless, there are broad («2000 km s" 1 wide), shallow absorption features at -23,000, -13,000, and -8000 km s _1 
in the C IV absorption region. Intervening systems at -79,100 and -75,500 km s _1 (offset from the system velocity) 
may be responsible for the "high-velocity" absorption, although we note in that case that the intervening absorber 
seems to have varied between epochs. The variation may also be due to intrinsic BAL components at high velocities. 

The broader absorption features in the C IV region at -13,000 and -8000 km s _1 are likely intrinsic. The higher- 
velocity component has weakened between epochs, while the lower-velocity component has disappeared altogether. 

1235+1453 

Our power law fit appears to match the continuum well in regions free of absorption and emission, so that the 
differences between epochs appear to be real changes in the absorption. The broad C IV absorber has deepened 
between epochs to extend out to -25,000 km s _1 . The absorber transmission is relatively constant from -13,000 out 
to -25,000 km s" 1 . 

1243+0121 

Similarly to 1235+1453, this source shows broad C IV absorption with a sharp onset at -5000 km s _1 . 

1314+0116 

This source is similar to 1235+1453 and 1243+0121 in that it shows deep, broad C IV absorption beginning at a 
detachment velocity of -4000 km s _1 . 

1331-0108 

C IV absorption from -15,000 to -23,000 km s _1 has apparently weakened gradually between epochs. 1331-0108 is 
a LoBAL, with broad Al III and Mg II absorption evident in the spectrum (not shown). 
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TABLE 1 
Observation Log 



LBQS 


SDSS 




z 


LBQS 


SDSS/HET 








B1950 


J2000 






Date 


Date 






(yr) 


0010 


- 0012 


001306.15 + 


000431.9 


2.16 


1987 Sep 25 


2000 Sep 6 


18.5 


-27.1 


4.1 


0018 


+ 0047 


002127.88 + 


010420.2 


1.82 


1988 Nov 2 


2000 Dec 7 


17.8 


-26.8 


1.3 


0021 


- 0213 






2.29 


1988 Sep 9 


2001 Dec 10 


18.7 




4.0 


0051 


-0019 


005355.15 - 


000309.3 


1.72 


1989 Sep 29 


2000 Aug 27 


18.7 


-27.1 


4.0 


0055 


+ 0025 


005824.75 + 


004113.3 


1.92 


1988 Sep 9 


2000 Aug 27 


17.5 


-28.6 


4.1 


0109 


-0128 


011227.60 - 


011221.7 


1.76 


1989 Sep 29 


2000 Sep 7 


18.3 


-27.3 


1.0 


1208 


+ 1535 


121125.49 + 


151851.5 


1.96 


1987 Mar 26 


2005 Apr 6 


17.9 


-27.4 


6.1 


1213 


+ 0922 


121539.66 + 


090607.4 


2.72 


1987 Mar 27 


2003 Feb 2 


18.1 


-28.0 


1.3 


1234 


+ 0122 


123724.51 + 


010615.4 


2.02 


1987 Apr 2 


2001 Feb 1 


18.0 


-27.7 


1.6 


1235 


+ 1453 


123736.42 + 


143640.2 


2.68 


1986 Dec 29 


2005 Mar 13 


18.6 


-27.7 


1.9 


1243 


+ 0121 


124551.44 + 


010505.0 


2.81 


1987 Apr 2 


2001 Jan 19 


18.5 


-28.3 


3.6 


1314 


+ 0116 


131714.21 + 


010013.0 


2.70 


1988 Mar 16 


2001 Mar 16 


18.6 


-28.0 


3.5 


1331 


- 0108 


133428.06 - 


012349.0 


1.88 


1988 Mar 21 


2002 Jun 1 


17.9 


-28.4 


1.9 



a The Bj magnitude given in the L BQS QSO catalog lIHewett et al.lll995l) . b The absolute i magnitude 
given in the SDSS DR5 QSO catalog dSchneider et alJI2007D . c The rest-frame time in years between LBQS 
and SDSS observations. 



TABLE 2 
Continuum and C IV BAL Data 



LBQS 


LBQS 


SDSS 


LBQS 


SDSS 


LBQS C IV 


SDSS/HET C IV 


Intervening System 


B1950 


L„(1400A) a 


L„(1400A) a 


L„(2500A) a 


L„(2500A) a 


EW (A) b 


EW (A) b 


Velocities 


0010 - 0012 


0.5 


1.4 


0.3 


1.6 


-7.3 + 1.3 


-3.8 + 0.8 




0018 + 0047 


1.4 


0.8 


1.2 


1.1 


-11.2 ± 1.2 


-13.5 + 0.5 


-27 (1.58), -57 (1.33) 


0021 - 0213 


1.8 




2.5 




-37.6 + 1.4 


-51.0 + 0.8 


0051 - 0019 


1.2 


0.5 


2.3 


1.5 


-40.6 + 1.1 


-29.1 + 0.8 


-29 (1.47), -56.2 (1.25) 


0055 + 0025 


4.7 


1.2 


5.1 


6.1 


-6.3 + 1.0 


-4.0 + 0.4 


-99 (1.07) 


0109 - 0128 


4.3 


1.2 


4.1 


1.9 


-12.6+ 1.0 


-29.5 + 0.6 


-60 (1.25) 


1208 + 1535 


1.5 


0.8 


1.1 


1.5 


-25.5 + 1.5 


-31.9 + 0.5 


-40 (1.59), -82.5 (1.23) 


1213 + 0922 


2.5 


2.3 


3.2 


3.8 


-13.8 + 2.0 


-9.5 + 0.5 


-42 (2.23) 


1234 + 0122 


1.3 


1.8 


1.6 


2.7 


-6.4 + 0.9 


-5.7 + 0.5 


-75.5 (1.33), -79.1 (1.30) 


1235 + 1453 


4.1 


1.0 


4.6 


1.7 


-29.2 ± 2.7 


-41.0 + 0.7 




1243 + 0121 


2.4 


3.0 


3.2 


1.7 


-34.4 + 2.5 


-34.0 + 0.5 




1314 + 0116 


2.0 


2.1 


1.5 


3.9 


-22.5 ± 1.9 


-30.9 + 0.5 




1331 - 0108 


0.9 


0.9 


2.3 


4.9 


-44.8+ 1.2 


-41.8 + 0.5 


-50.7 (1.43) 


a Monochromatic luminosities 


are given in units of 10 31 erg 


s — 1 Hz^ 1 and 


include a factor of (1 + z) "bandpass 


correction" (Hoeej|1999|). 



0021-0213 was not observed with the SDSS, and we do not include it in our continuum variation studies.* 5 In the C IV region, 1 A is 
ss200 km s _1 . c Velocities of any intervening systems are given in units of 1000 km s 1 offset from the rest frame. The corresponding rcdshift 
for each system is given in parentheses. 
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TABLE 3 
CIV Balnicity Indices 51 



t t~) /^i c T~> 1 ncn 

LBC^o BiyoU 


HF03 


LBQo Bl 


bubb/ Ha 1 Bl 


LBCjb BlQ 


bDbb/ Hh 1 B1q 




BAL? 


(km s *■) 


(km s ) 


(km s ) 


(km s *■) 


001 — D01 2 

UU1U uUli 


No 


18 





48 


o 


nm s -i- nn47 

uuio uui j 


Yes 


10 


39 


309 


257 


0021 - 0213 


Yes 


4042 


6108 


4042 


6108 


0051 - 0019 


Yes 


2474 


191 


2474 


191 


0055 + 0025 


No 


61 





61 





0109 - 0128 


Yes 


746 


1939 


746 


1939 


1208 + 1535 


Yes 


2517 


3302 


2517 


3302 


1213 + 0922 


No 


47 


25 


47 


25 


1234 + 0122 


Yes 














1235 + 1453 


Yes 


1342 


4820 


1342 


4820 


1243 + 0121 


Yes 


3622 


3764 


3622 


3764 


1314 + 0116 


Yes 


257 


3077 


257 


3077 


1331 - 0108 


Yes 


6568 


5977 


6568 


6199 



a Bl and BIq were calculated for the LBQS, SDSS, and HET spectra using our continuum 
fits as described in fl2.2| 



TABLE 4 
Varying Regions 



LBQS 




Av 


(t»> 


B1950 


(km s ) (km s ) 


(km s 


0018 + 0047 


2300 


2300 


-3600 


0021 — 0213 


6000 


2500 


—24000 






3500 


— 16700 


0051 - 0019 


9800 


5000 


-20700 






3000 


-12200 






1800 


-8400 


0055 + 0025 


2400 


2400 


-12700 


0109 - 0128 


10500 


1700 


-23000 






2500 


-17600 






6400 


-8600 


1208 + 1535 


10800 


3700 


-20300 






1800 


-17300 






2600 


-14100 






2600 


-10600 


1213 + 0922 


1200 


1200 


-21900 


1234 + 0122 


2000 


2000 


-7900 


1235 + 1453 


9700 


2500 


-26700 






1700 


-22800 






1700 


-20700 






1400 


-18900 






1200 


-12300 






1200 


-10500 


1243 + 0121 


2600 


1200 


-16400 






1400 


-12800 


1314 + 0116 


6100 


1700 


-21100 






1200 


-15400 






1600 


-11300 






1600 


-6800 


1331 - 0108 


3500 


2300 


-21200 






1200 


-19300 



a XIa-i; is the sum of velocity widths of all varying 
regions in a source, so only one value is given per 
source. For sources not listed in this table, we found 
no varying regions, and Sat; — 0. 
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TABLE 5 
Trough Onset 
Acceleration Limits 



LBQS 


Upper Limit 


B1950 


(cm s — 2 ) 


0018 + 0047 


0.14 


0021 - 0213 


0.15 


0051 - 0019 


0.15 


1235 + 1453 


0.12 


1243 + 0121 


0.17 


1314 + 0116 


0.17 


1331 - 0108 


0.12 



. Lundgren 
et al. (2007) 



o 
o 
o 



Barlow (1993) 



0.1 



At ( 



(years ) 



This 
Study 



1 



Fig. 1. — A comparison of the rest-frame EWs (averaged over two epochs) of the C IV BALs for the sources in B93 (triangles), L07 
(crosses), and our study (squares). For each source, —(EW) is plotted against the rest-frame time between epochs, A.t sys . Vertical bars 
mark the error in EW. Note that the x-axis is logarithmic. 





LBQS 0010-0012 - 
1987 Sep 25 ; 




SDSS J001306. 15+000431. 8: 
2000 Sep 6 : 



1500 2000 2500 



Fig. 2. — Plot of the LBQS (top) and SDSS (bottom) spectra of each source. The best model fit for the continuum and strong emission 
lines in each epoch is plotted with a thick red line. The data are binned by a factor of of 2 over the instrumental sampling. The y-axis is 
the observed-frame flux density in units of 10 — 17 erg s _1 cm -2 A -1 at the rest-frame wavelength given by the ir-axis. In cases where the 
continuum in the C I V absorption region is not well-constrained, we have used the LBQS fit (with data extending to shorter wavelengths) 
as a guide. See £|2.1I for additional details. The sources shown and their dates of observation are indicated in the plots. Figures for the 
remaining sources are available in the electronic edition of the Journal. The printed edition contains only a sample. 
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0010-00,12 




l v ' I -I 1 I yMrt i ^-^J. 



0055+0025 




"I ' ' i t 

1208+1535 



-3x10 4 -2.5x10 4 -2x10* -1.5x10 4 - 1 4 

Velocity (km s 1 ) 



-5000 



1213 + 0922 





1331-0108 



-3x10 4 



-2.5x10 4 




-2x1 4 



-1.5x10 
Velocity (km s 1 ) 



-5000 



Fig. 3.— Ratio spectra R(v) (as defined in fl2~T) of C IV BAL regions for the BAL QSOs in our sample. The spectra have been divided by 
the best-fit reddened power law with C IV and Si IV emission lines for each epoch. LBQS spectra are shown as thin black lines (with error 
bars), while the SDSS and HET spectra are shown as thick lines. The SDSS spectra have been convolved with a Gaussian to approximate 
the LBQS resolution. Vertical tick marks above the spectra indicate locations where lines from intervening systems may be present (see 
§??). The thick, horizontal lines above the spectra mark regions at least 1200 km s — 1 wide which varied between epochs, as described in 
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Fig. 4. — Top panel: The horizontal bars indicate velocity regions which were observed to vary (in the sense of 93.1.10 for the sources 
listed. Bottom panel: The number of times a particular velocity was included in a varying region. 



1000 



2000 



3000 



4000 



5000 



Av (km s ) 

Fig. 5. — The number of regions in all BAL QSO sources with width A?; that varied between epochs. We only consider varying regions 
with widths >1200 km s -1 wide in order to minimize false detections due to statistical variation. We do not consider varying regions of 
smaller widths. 
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0.1 



0.3 



0.4 



0.5 



<AR> 

Fig. 6. — The distribution of (AiJ}|, the average change in the ratio spectrum between epochs in regions at least 1200 km s wide 
which varied (as defined in £|3.1, It . Our definition of "varying region" is not sensitive to very shallow absorption variation, so there may 
be additional variation with small | ( AR) | that we have not identified. 




Fig. 7. — Top panel: The distribution of (R(v)) in the more deeply-absorbed epoch. Middle panel: The distribution of (R(v)) in the 
less-absorbed epoch. Bottom panel: The distribution of (R(v)) in all 1 A bins from 1400-1550 A for both epochs of all sources in our 
sample. The normalization of the bottom panel is arbitrary. See ^3.1.31 for more information. 
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0.5 1 



Fraction of Bins With R(v) S 0.8 

Fig. 8. — The fraction of 1 A bins in the C IV absorption region with ratio spectra varying by Ai?(A)j > 0.2 is correlated with the 
fraction of BAL bins (R(X) < 0.8) in that region. The best linear fit (Equation [TJ is shown. 




Fig. 9.— AEW vs. (EW) (in km s" 1 ) across the entire C IV absorption region for each BAL in our sample. The symbols described in 
the legend indicate the samples of B93, L07, and this study. 
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F„(2500 A) 



F (1400 . 



-0.5 0.5 

AJog10(F v ) 

Fig. 10. — The change in C IV region EW as a function of the change in flux. The error on the flux is assumed to be dominated by the 
LBQS flux error, estimated to be about 15% (Hcwctt ct al. 1991). The EW errors are calculated from the formal error on EW, assuming 
the continuum error is 5%. Top panel: the change in flux, Alog(F l ,), is calculated at 2500 A. Middle panel: Alog(.F„) is calculated at 
1400 A. Bottom panel: Alog(.F„) at 1400 A (y-axis) is plotted against Alog(i ;, J/ ) at 2500 A (rr-axis). The two are strongly correlated. 




At (years) 

Fig. 11. — The change in EW against the rest-frame time between epochs. The envelope of AEW clearly expands with time. We show 
data from B93, L07, and our study, identified with symbols as shown in the legend. The x-axis is logarithmic. The solid lines indicate the 
square root of the unbiased sample variance calculated from a sliding window of 15 time-ordered entries. 
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At (years) 

Fig. 12. — The magnitude of fractional change in EW against the rest-frame time between epochs for sources in our study, as well as for 
those of B93 and L07. We have constructed a matched sample of sources by including only the B93 and L07 sources with j (EW) | greater 
than the minimum in our sample. The x-axis is logarithmic. The symbols described in the legend indicate the samples of B93, L07, and 
our sample. 
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